EFFECTIVE STRING PICTURE FOR CONFINING GAUGE 
THEORIES AT FINITE TEMPERATURE 



M. CASELLE AND M. PANERO * 

Universita di Torino and INFN, Sezione di Torino, 
Via P. Giuria, 1 
1-10125 Torino, Italy 
E-mail: caselle@to.mfn.it, panero@to.infn.it 

M. HASENBUSCH 

NIC/DESY Zeuthen 
Platanenallee 6, 
D-15738 Zeuthen, Germany 
E-mail: Martin.Hasenbusch@desy. de 



We study the "effective string picture" of confinement, deriving theoretical predic- 
tions for the interquark potential at finite temperature. At low temperatures, the 
leading string correction to the linear confining potential between a heavy quark- 
antiquark pair is the "Liischer term". Assuming a Nambu-Goto effective string 
action, subleading contributions can be worked out in an analytical way. We also 
discuss the contribution given by a possible "boundary term" in the effective ac- 
tion, and compare these predictions with high precision results from simulations 
of lattice Z2 gauge theory in three dimensions, obtained with an algorithm that 
exploits the duality of the Z2 gauge model with the Ising spin model. 



1. Effective String Picture of Confinement 

The string picture is an effective framework which is expected to provide 
a good physical description of the infrared behavior of confining gauge 
theories. The basic idea is simple: two confined color charges behave as if 
they were joined by a thin flux tube, which can fluctuate like a vibrating 
string. The dynamics of the world sheet spanned by the string during its 
time-like evolution is described by an effective action S o{{ , and for massless 
string fluctuations, the simplest choice for a candidate effective action is 
the Nambu-Goto string action: S clt = a ■ A, which is proportional to the 
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area A of the world-sheet surface; a is the string tension, appearing as the 
parameter of the effective theory. 

For a three-dimensional system with extension L 2 xL (with L s ^> L) and 
periodic boundary conditions along the short direction, the temperature T 
is proportional to l/L, and in that case the result for the expectation value 
of the Polyakov loop correlation function reads: 

-aRL+k 

where rj is Dcdekind's function. The term associated with the minimal 
world sheet surface induces the exponential area-law falloff responsible for 
the linear rise in the interquark potential V(R), while the first non-trivial 
contribution in S ot{ results in the determinant of the Laplace operator, and 
the corresponding contribution to the interquark potential V(R) — in a 
regime of distances shorter than ^ — is the Luscher term 1 : 

V(R) = -i ln(pt (jR )P( )) ~ aR - ^ (2) 

Inclusion of further terms 2 in the expansion of the world sheet area 
results' 1 in a contribution involving a combination of Eisenstein functions 3 : 
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However, such a contribution is still under debate 2 . 

On the other hand, it is also possible to include a "boundary term b " in 
the effective action: a perturbative expansion in b (a parameter proportional 
to the coefficient of the "boundary term" in the effective action), induces a 
leading order correction like: 

with a short distance contribution to V(R) reading: — -^p- 



a This calculation involves a C, function regularization. 

b Such a boundary term is related to derivatives of the h field (which describes transverse 
displacements with respect to the minimal area surface of the world sheet), evaluated 
along the Polyakov lines 
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2. The Model: Z 2 Lattice Gauge Theory 

We run numerical simulations of the Z 2 lattice gauge theory in three space- 
time dimensions. This choice has various motivations: the effective string 
picture is believed to be independent of the underlying gauge group; the 
Z 2 gauge group is interesting from the perspective of the center role in 
confinement ; the reduced configuration space of this theory and its duality 
with respect to the Ising spin model enable one to get high precision results 
within a reasonable amount of CPU time. 

The pure 3D lattice gauge model is described in terms of a x ^ variables 
(taking values in Z2) defined on the lattice bonds; the dynamics is governed 
by the standard Wilson action, which enjoys Z 2 gauge invariance d . The 
partition function reads: 



and the system may exist in different phases: a confined, strong coupling 
phase, with massive string fluctuations for (3 < 0.47542(1) 4 ; a confined, 
rough phase, with massless string fluctuations (this is the regime we studied 
in our simulations); a deconfined phase for j3 > 0.7614134(2) 5 . 

This model is dual with respect to the Z 2 spin model in 3D, and we 
exploited this property to express a ratio between Polyakov loop correlators 
G{R) = (Ft (-R)P(O)} of the gauge theory as a product of expectation values 
of one-link variables in the modified spin ensembles. A similar algorithm is 
the so-called snake algorithm 6 . 

We used multi-level updating and a hierarchical organization of sub- 
lattices, and the CPU time turns out to be roughly proportional to the 
inverse temperature L, and virtually independent of the distance R between 
the quark sources, thus the algorithm is particularly useful in a regime of 
very large interquark distances. 

3. Numerical Results 

Let F(R, L) be the free energy associated with the presence of a heavy 
quark-antiquark pair at finite temperature: G(R) = e~ F ( R,L \ We stud- 
ied "quantum contributions" in free energy differences, by measuring the 



c 7i2 is the center of continuous gauge groups like SU(2) or Sp(N). 

d 7i2 gauge transformations act as local flips of c Xtll variables living on the lattice bonds 
which meet at a given site. 
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following quantity: 

Q(R, L) = F(R + 1, L) - F(R, L) - aL 



(6) 



T=2T_c/3 



a 




Figure 1. Q(R, L) for L = 12 (i.e. T = 2T c /3) at /3 = 0.75180. z is defined as: z=^. 
Solid curves correspond to the free bosonic string prediction, while dashed lines include 
the LO Nambu— Goto correction. Pure area— law corresponds to Q = 0. 



Fig. 1 shows that at "high temperatures" (L < 2R) our numerical 
results are in good agreement with the NLO prediction from the Nambu- 
Goto string, while a pure area law is definitely ruled out, and the LO term 
alone is not sufficient to describe the data. We also found that the coefficient 
of a possible "boundary term" for this model seems to be very small, likely 
zero. 

As it concerns the (L > 2R) regime, Fig. 2 shows the deviation 2 of a 
quantity proportional to [Q{R, L) — Q(R — 1 , L)] from the free string pre- 
diction. The chosen normalization allows a meaningful comparison among 
different LGT's in 3d: SU(2) gauge theory 7 (crosses), SU(3) gauge theory 8 
(white squares), and Z 2 gauge theory 2 (black squares). The three models 
display the same qualitative behavior, and, in particular, the data for Z 2 
and SU(2) (which are groups with the same center — namely: Z 2 itself) 
are compatible within errorbars. This may be a signature of the relevance 
of center degrees of freedom to the confinement mechanism. 
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Comparison of SU(2), SU(3) and Ising 
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Figure 2. Behavior of different gauge models in the low temperature regime. 

4. Conclusions 

We studied confining gauge theories at finite temperature, and tested the 
theoretical predictions of the Nambu-Goto effective string for Z 2 lattice 
gauge theory, both at large and short interquark distances. Our algorithm 
exploits the duality of the model, and this enabled us to explore a wide 
range of distances, detecting next-to-leading order effects. Our data seem 
to rule out a "boundary term" in the effective string action describing 
the present gauge model. Finally, we also made a comparison with some 
different gauge models. 



References 

1. M. Liischcr, K. Symanzik and P. Wcisz, Nucl.Phys. B173 (1980) 365. 

2. M. Caselle, M. Hasenbusch and M. Panero, JHEP 0301 (2003) 057. M. Caselle, 
M. Panero and P. Provero, JHEP 0206, 061 (2002) and references therein. M. 
Caselle, M. Hasenbusch and M. Panero, in preparation. 

3. J. P. Serre, "A course in Arithmetic", Springer- Verlag, New York, 1980. 

4. M. Hasenbusch and K. Pinn, J.Phys. A30 (1997) 63. 

5. H.W.J. Blote, L.N. Shchur, A.L. Talapov, Int.J.Mod.Phys. CIO (1999) 1137. 

6. Ph. de Forcrand, M. D'Elia and M. Pepe, Phys. Rev. Lett. 86 (2001) 1438. Ph. 
de Forcrand, M. D'Elia and M. Pepe, Nucl.Phys. Proc.Suppl. 94 (2001) 494. 

7. M. Caselle, M. Pepe and A. Rago, in preparation. 

8. M. Liischer and P. Weisz, JHEP 0207 (2002) 049. 



